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1. The Yang-Mills gradient flow 

The basic idea of the Yang-Mills flow consists in introducing an extra dimension to our gauge 
fields, that we will call flow time and denote by t (with dimensions [length^]), and study how fields 
evolve according to the flow equation 

dBii(x,t) 

^ =DyGy^{xfl), ( 1 . 1 ) 

with initial condition B^{x,t = 0) = A^{x), and G^y being the field strength 

Gv^{x,t) = dvB^{x,t)-dvB^{x,t) + [Bv{x,t),Bfi{x,t)]. (1.2) 

Although it has been used in different contexts in our field [1, 2, 3], the recent growth of 
works and interest in the gradient flow have to do mainly with two facts: first the well understood 
properties of composite operators defined at positive flow time (t > 0) under renormalization and 
its continuum limit [4, 5, 6], and second with the concrete proposals of using the gradient flow 
to give a precise definition of the renormalized coupling [4], and as a general tool to renormalize 
composite operators [5, 6]. In this contribution I will try to summarize the recent works of the 
community in these directions. 

The key idea behind most of these applications is that the r.h.s. of the flow equation (1.1) is 
nothing but the gradient of the Yang-Mills action, and therefore the gauge field becomes smoother 
along the flow: the flow smears the gauge field over a region of radius killing the UV diver¬ 
gences and making composite gauge invariant observables automatically renormalized. 

For example the energy density, an observable that we will see throughout all this proceedings 
contribution, defined as 

{E{x,t)) = -^{TrGf,v{x,t)G^v{x,t)), (1.3) 

is finite for t > 0. This is not an accident of this particular observable, or of leading order perturba¬ 
tion theory. The existence of a continuum limit (to all orders in perturbation theory) for all gauge 
invariant composite operators at positive flow time was proved in [5]. 

Since t^{E{t)) is a dimension-less renormalized quantity that depends on a scale (given by 
\/^), it is a natural candidate for scale setting and to define the renormalized coupling, as we will 
see in the following pages. 

The second application that we will cover is related with the behavior of an operator like 
{E{t)) as t —)• 0. The correlation function will become singular in a way that can be described by 
an expansion in renormalized operators with singular coefficients [5] 

E{x,t) = ^c«(0{(9“}«W + m , (1-4) 

a 

where { 0 ^} r { x ) are the set of renormalized operators that “mix” with E{t), and Ca{t) are the 
singular coefficients. This so called small flow time expansion, has recently been used in the context 
of thermodynamics, as has been reviewed in this conference [7]. But the idea of the small flow time 
expansion has a more general scope, and we will dedicate some time to this matter. 

For any of these applications taking the continuum limit of flow quantities is a fundamental 
step. The community has recently dedicated some effort in understanding and controlling the cutoff 
effects of flow quantities, and these will also be reviewed here. 
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2. The gradient flow coupling 


For the moment, most of the applications of the gradient flow to the definition of a renormal¬ 
ized coupling make use of the energy (or action) density at positive flow time 

{E{x,t)) = -^{TvG^^ix,t)G^y{x,t)). (2.1) 

As we have said in the previous section, {E{x,t)) is a renormalized quantity, and in fact it has been 
shown in [2] that ^(.r,?) has a perturbative expansion given by 


where ci is a finite coefficient. 

This immediately suggests two applications [2]. First one can define a scale to fixed by the con¬ 
dition tQ{E{x,tQ)) = 0.3. This subject was reviewed extensively in the last lattice conference [8], 
and we will not talk more about scale setting. Second, one can use t^{E{x,t)) to define a renormal¬ 
ized coupling 




1671 ^ 


t^{E{x,t)) 


f=l/8/j2 ■ 


(2.3) 


Although Eq. (2.3) can be used as it is, it requires the existence of a window l/L<Cft<Cl/a 
so that the determination of is free of both lattice artifacts and finite volume effects. This 

Windowing problem can be solved in a very elegant way by relating the renormalization scale with 
the finite size of the system (finite-size scaling) [9] 


At = 


1 


(2.4) 


where c is a constant that relates both scales. Usually, in an abuse of notation, one writes for 

this running coupling. Since this coupling naturally feels the finite volume in which it is defined, 
the coefficients in the perturbative expansion Eq. (2.2) will depend on the choice of boundary 
conditions, as will the concrete coupling definition. All these coupling definitions can be written as 


glp{L) = ^ h'^{E{x,t)) 


y/^—cL 


(2.5) 


where JY will depend on the boundary conditions and choice of c. 

Running coupling schemes have been proposed with periodic [10], Schrddinger functional 
(SE) [11], twisted (a la t’Hooft) [12], and mixed open-SE [13] boundary conditions. The advan¬ 
tages of each particular choice for a concrete project have to be evaluated taking into account the 
following points 


Cutoff effects Schemes that break the invariance under translations, like the SE or open-SE schemes 
require further work to ensure an scaling towards the continuum. Quantum field theories 
defined in manifolds with boundaries require additional boundary counterterms for renor¬ 
malization [14]. Eor the concrete case of Yang-Mills theories on the lattice this implies that 
additional boundary counterterms have to be introduced to avoid G{a) cutoff effects [15]. 
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Moreover, Wilson fermions require additional boundary and bulk counterterms to ensure 
^{a^) scaling [16, 17]. 

On the other hand schemes that preserve the invariance under translations, like the ones that 
use periodic or twisted boundary conditions, do not have these extra difficulties. Automatic 
ff{a) improvement is guaranteed when working with massless Wilson quarks [18, 17] with¬ 
out having to introduce any counterterms. 

Zero modes in perturbation theory One important characteristic of the schemes with either twisted, 
SF or open-SF boundary conditions is that the field configuration that minimizes the action 
is unique up to gauge transformations. On the other hand this is not the case for the periodic 
case. The presence of zero-modes that are not gauge degrees of freedom complicates substan¬ 
tially any perturbative computation that one wants to address [19]. Moreover the observable 
t'^{E{x,t)) has a non-analytic expansion in otj,^, leading to a non-universal j3—function. 

At this point is important to make clear that this is a problem of perturbation theory, and 
any non-perturbative computation is insensitive to this issue. But one should not under¬ 
estimate the usefulness of perturbative computations in order to get insight into the full 
non-perturbative results. Also, we recall that many applications of these running coupling 
schemes need at some point to match with pertubation theory. This matching is not only 
more difficult because of the additional complications of the perturbative computations, but 
one might have to do it at higher energy scales than in the case of schemes with a universal 
j8—function. 


Critical slowing down Recently the well known problem of topology freezing [20,21] that affects 
large volume simulations has been found to affect also studies of the running coupling [22]. 
Although in [22] a way to overcome this problem is proposed, essentially by restricting 
the measurements to a sector of fixed fopological charge, the open-SF scheme is the only 
scheme that address this issue in a more simple and elegant way by choosing some boundary 
conditions that allow the topological charge to flow in and out of the lattice [13]. 


In all these studies a central role is played by the step-scaling function 


<Ts{u) — gGFi^^) 


glp{L)=u ’ 


( 2 . 6 ) 


that measures the change of the coupling when the renormalization scale is changed by a factor 
\/s. One should note that this step-scaling function can be computed very easily on the lattice by 
changing the lattice size Lja while keeping the bare parameters fixed. 


2.1 Comparison of the gradient flow coupling and the SF coupling 

Studies of the running coupling using the idea of finite-size scaling have been done until re¬ 
cently using the so called SF coupling [15]. In this case one imposes SF boundary conditions on 
the fields and measures how the action depends on the induced background field. Here we will 
review the differences between the SF and the GF couplings. 

In Tab. 1 we can see some data [11] of both the SF and the GF couplings. These data corre¬ 
spond to 77/ = 2 simulations at constant volume L\ ~ 0.4fm, and therefore at constant renormal¬ 
ization scale jX. This is why the values of ggp{Li ) are constant within errors. On the other hand the 
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L/a 

6 

8 

10 

12 

16 


5.2638 

5.4689 

5.6190 

5.7580 

5.9631 

^sea 

0.135985 

0.136700 

0.136785 

0.136623 

0.136422 

^meas 

12160 

8320 

8192 

8280 

8460 


4.423(75) 

4.473(83) 

4.49(10) 

4.501(91) 

4.40(10) 

?gf(m)(c = 0.3) 

4.8178(46) 

4.7278(46) 

4.6269(47) 

4.5176(47) 

4.4410(53) 

^gf(m)(c = 0.4) 

6.0090(86) 

5.6985(86) 

5.5976(97) 

5.4837(97) 

5.410(12) 

^gf(m)(c = 0.5) 

7.106(14) 

6.817(15) 

6.761(19) 

6.658(19) 

6.602(24) 


Table 1: Comparison between the gradient flow coupling and the Schrbdinger functional cou¬ 
pling g|p in a set of Nf = 2 ensembles at constant volume given by Li = 0.4fm [11]. 

values of the GF coupling are not constant within errors. This is a manifestation of the cutoff effects 
in the GF coupling, since the values of gQp{L) of some particular ensemble depend not only on L, 
but also on the lattice spacing a. We will talk about cutoff effects in more detail later, but for the 
moment it is enough to say that the relative size of cutoff effects decreases when c increases [11]. 

Tab. 1 also shows that the GF coupling is much more precise than the SF coupling. We have to 
take into account that the numbers of the GF coupling were produced with one order of magnitude 
less statistics than the ones with the SF coupling. In this sense the observable used to define the GF 
coupling has much smaller variance, and therefore leads to more precise determinations. The data 
in Tab. 1 also show that this variance increases with the value of c. A more quantitative comparison 
of the variance of both coupling definitions can be found in [23] were it is argued that the variance 
of the GF coupling is around 50-100 times smaller at these volumes. 

An important point is how this variance scales both with the lattice spacing a and the physical 
volume L. Summarizing the studies made in [24, 23], we can say that the variance of the SF 
coupling diverges as one approaches the continuum like 1/a, while the variance of the GF coupling 
is roughly independent of a in the scaling region. On the other hand, and as R. Sommer realized 
very soon, the variance of the SF coupling behaves in small volumes (small values of the coupling) 
as ggp, while the GF coupling behaves as gQp. These different behaviors are used by the ALPHA 
collaboration to devise an optimal strategy to compute the running coupling in QCD [25]. In this 
approach the SF coupling is used in the small volumes to match with perturbation theory, and the 
GF coupling is used in the larger volumes to match with a hadronic scheme. 

3. Cutoff effects and the gradient flow 

On the lattice the gradient flow is defined through the equation 

a^dtVi,{xd) = -gl{r'dl^S^{V)}V^{xd ), V^{x,f)) = U^{x ), (3.1) 

where (x) are the links associated with the fundamental gauge field, {x,t) are those associated 
with the flow field, and 5'fl(F) is any lattice action (not necessarily the same used in the simu¬ 
lation). The derivative has to be understood as Lie-algebra valued (for more details see for 
example [4]). Eq. (3.1) is simply a lattice version of the flow equation Eq. (1.1), in which the links 
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are evolved according to the gradient of a lattice action. It has become clear recently [26, 27], that 
the choice of the action that defines the flow, together with the choice of action for the simulation of 
the fundamental gauge field and the choice of discretization for the observable are the ingredients 
that determine the cutoff effects on any flow quantity. 

Step-scaling studies and computations of the running coupling have only one systematic effect 
that needs to be taken under control: taking the continuum limit. In the studies with the SF coupling 
it is not strange to use small lattices, including those with L/a = 4,6. On the other hand schemes 
that use the gradient flow coupling typically use larger lattices. This difference does not mean 
that gradient flow observables have large cutoff effects per se, since the natural scale that controls 
the size of cutoff effects for flow quantities is not the whole lattice (ajL), but the region that is 
smeared by the flow {a/s/^ = ajcL). This suggests that larger values of c might lead to smaller 
cutoff effects, and this is in fact the case [11]. But this comes at a price: first, as we have seen, 
couplings defined at larger values of c have larger errors. In addition, in schemes where the lattice 
has boundaries in the time direction (SF or open-SF), the boundary effects become more noticeable 
at larger c [13]^ 

In order to illustrate the dangers of the continuum extrapolation in step scaling studies, I will 
use some data from these proceedings contribution [28] (see Fig. 1). As we can see, when two 
different discretizations for the observable E{x,t) are used, the continuum extrapolations of the 
step scaling function do not agree for small values of the parameter c (Fig. 1 (a)). Note that in this 
particular case, if one only looks at one observable and only at one value of c, it would be difficult 
to spot the problem. In particular the continuum extrapolations look reasonable even for c = 0.25. 
On the other hand for large enough c ~ 0.4, both continuum extrapolations agree nicely (Fig. 1 
(c)). Figs. 1 (b) and (d) show that the disagreement is larger at strong couplings. 

This example illustrates very well the reasons behind the efforts that different groups have put 
into understanding cutoff effects of flow quantities. 

3.1 Tree-level correction 

There are many possible choices for lattice gauge actions, but a general class of them can be 
written symbolically as 



These include the Liischer-Weisz tree-level improved action (co = 5/3,ci = — 1 / 12 ,C 2 = C 3 = 
0) and the Iwasaki action (cq = 3.648,ci = —0.331,C 2 = C 3 = 0). 

One has to choose several actions to study a flow quantity, first the action used for the simula¬ 
tion (i.e. to produce the configurations), second the action used in the flow equation (3.1). Finally, 
note that {E{x,t)) is nothing more than an action density, and therefore, if interested in this observ¬ 
able, any particular discretization of {E{x,t)) can be understood as a third choice of action. Each 
of these actions can be different, and therefore we have several sets of coefficients a: one for the 
action, one for the flow, and one for the observable. 

'Note that this can be ameliorated by simply increasing the ratio T/L. 
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Open symbol: plaquette; Filled sumbol: clover 



(8 -» 16, 10 ^ 20, 12 24) lineai' continuum extrapolation 



(a) 


(b) 


Open symbol: plaquette; Filled sumbol: clover (8 16, 10 -> 20, 12 —* 24) linear continuum extrapolation 



(C) (d) 

Figure 1: Continuum extrapolation of the step scaling function G 2 {u) and running coupling us¬ 
ing different discretizations for the observable E{x,t). Top row: c = 0.25, Bottom row c = 0.4. 
Left Column: Continuum extrapolation for 02(5). Right column: 02 {u)/u in the continuum as a 
function of u. See the text for more details. (Source [28]). 


The influence of these coefficients to tree-level has been clarified recently [26, 27]. The per¬ 
turbative expansion of our observable of interest reads 

{E{x,t)) = ^oit,alL)g^ + ^(/), (3.3) 


where Sq can be computed on the lattice. In fact in the case of an infinite lattice we have 


t^S’Q{t,alL) = 


128712 


1 + ^ 
t 


j(«) _ _ 3 j(/) 


+ 0’{a^ 


(3.4) 


where the coefficients parametrize the tree-level cutoff effects produced by the observable 

discretization, the action used in the simulation and the action used to define the flow respectively. 
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The values of these parameters for some popular choices of discretizations are 


dio.aJ) 


' -1/24 
< 1/72 
_ -5/24 


plaquette 
Liischer-Weisz . 
Clover 


At this point there are some important observations to make 


(3.5) 


1. The particularly popular choice of using the clover discretization to define flow quantities 
(couplings, or the scales Iq^wq, ...), does not seem to be a very good idea from the perspec¬ 
tive of perturbation theory. This discretization is the one with larger tree-level cutoff effects: 
5 times larger than the plaquette definition, and 15 times larger than the Liischer-Weisz dis¬ 
cretization. 

2. There is a cancellation between cutoff effects produced by the observable and those produced 
by the action and the flow. For example the choice of a Wilson action, Wilson flow, and 
clover observable leads to rather small tree-level cutoff effects. This is due to an accidental 
cancellation of the cutoff effects produced by the action and flow on one hand, and those of 
produced by the observable on the other. Note that this accidental cancellation disappears 
as soon as the setup is changed. For example, using the Symanzik flow and Liischer-Weisz 
action, but keeping the clover observable leads to a very poor scenario, with tree-level cutoff 
effects around 7 times larger. 


The perturbative computation of Eq. (3.4) immediately suggests an improved definition for the 
coupling 


glp(L) = ^ h^(E(x,t)} 


y/^—cL 


(3.6) 


where ^ (i.e. Using the tree-level lattice computation instead of the continuum result to 

define the coupling). This definition of the coupling has no cutoff effects to tree-level, and has been 
used in most step scaling studies (see [11, 12, 13, 28, 29]). The exception were the studies using 
periodic boundary conditions (like [10, 30, 31]), or smeared gauge actions (for example [34]), since 
the perturbative computation in these setups is more involved (see the recent work [26] for the case 
of periodic boundary conditions). 

This improved definition of the coupling has a big impact on the cutoff effects. A nice example 
is given in [26], where the step-scaling study of [10] is repeated with the improved definition of the 
coupling Eq. (3.6). The result can be seen in Eig. 2. This example is specially impressive, partially 
because the original data of [10] has a choice of discretizations (Eiischer-Weisz action, Symanzik 
flow and Clover observable) that, as we have commented, make the tree-level correction specially 
large. 

The work [10] also proposes to use a general class of actions given by Eq. (3.2) setting C 2 = 
C 3 = 0. Since the two remaining coefficients are related by cq -|- 8ci = 1, each choice of action 
has one free parameter, resulting in three parameters that can be chosen arbitrarily. The concrete 
proposal is to choose them to cancel the tree-level cutoff effects of {E{x,t)) in infinite volume to 
order a^,a‘^ and a^. It is the opinion of the author that this criteria to choose the coefficients is not 
very appealing, especially taking into account that one can remove all the tree-level cutoff effects 
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Nf=4 s=1.5 step function tree-level improved 



Figure 2: Comparison of the continuum extrapolation of the step scaling function of the coupling 
defined with the continuum norm and the tree-level improved definition. (Source [26]). 


(to all orders in a) with the simple coupling definition of Eq.(3.6), and there is no good reason to 
think that this choice of coefficients will result in reduced cutoff effects for any other quantity or to 
higher orders in go- 


3.2 Phenomenological improvement 

In any 0’{a) improved setup, the lattice gradient flow coupling can be expanded as 

VHt=cL = 


gQp{L-,a/L) = ^ h^{E{x,t)) 


where, as we have seen in fhe previous section, A{c,L) depends nof only on fhe volume and fhe 
parameter c, buf also on fhe choice of boundary conditions, fhe acfion, fhe flow, efc... The proposal 
of [30] is fo subfracf fhe cutoff effecfs a^A{c,L) by modifying fhe definition of fhe coupling to 


gQp{L-,ajL) = ^ ^t^{E{x,t-{-a^r)) 


\/^—cL 


(3.8) 


The paramefer T is chosen fo cancel fhe leading cutoff effecfs. One can see fhis r-shiff as adding an 
term to the coupling definition 


g^^{L\a/L) = gQp(L;a/L) + -j^{t'^dtE{xd)) 


y/^—cL 


+ ^(a'^), 


(3.9) 


Recently a similar approach in spirit has been taken in [32]. In this case instead of canceling 
the term by a shift in the flow lime, fhe aulhors use a linear combinafion of fhe coupling 

defined wilh differenl discrelizalions of fhe observable 


gQ^{L-,a/L) — XgQp.QQyg^(L;a/L) + (1 — 2f)gQp.piaquette(^i‘^/^) j (3.10) 
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2 

0 1/18M/16- 1/12^ 


c=0.2, s=2 


ro = 0.0 I— ▼— 

TO = 0.02 I— ■- 

To = 0.04 I— #- 
m = 0.06 >— 


(a/Lf 


OC 

I 


(a) 


1.2 

1 

0.8 

0.6 

0.4 

0.2 

0 

0 1/20^ 1/16^ 1/12- 


g2(L)=1.8,ro = -0.02 

g-{L) = 2.2, TO = -0.02 . 
g-{L) = 2.6, To = -0.02 I 
gHi-) = 3.0,ro = -0.01 . 
g\L) = 2.4,To = -Om I 


(a/Lf 

(b) 


Figure 3: Left: Continuum extrapolation of the step-scaling function for different values of T 


where the parameter X is chosen to eliminate the leading cutoff effects. 

Obviously these are always legitimate things to do, that cannot change the values obtained af¬ 
ter a properly done continuum extrapolation. In this context one should avoid giving a dependence 
to T or X on the bare parameters (like j3), since this would result in a different shift for each 
value of a, and could give the wrong continuum extrapolations. One should also be careful with 
extrapolating to the continuum with a constant fit, even if the data “looks fiat” after having applied 
the correction, since one could artificially reduce fhe statistical error of the extrapolation. Fig. 3 
shows how different values of T results in data with different cutoff effects, but whose continuum 
extrapolations are in agreement. Moreover these recipes are simple to implement since the inte¬ 
gration of the flow equafions naturally gives {E{x,t)) at many values of t, and the evaluation of a 
different observable is always cheap. 

In the opinion of the author the weak point of these approaches is that the improvement of 
flow quantities requires one to improve the flow equation, the discretization of the observable and 
the action used to generate the configurations [8, 27]. A single parameter (t or X) is not able to 
achieve this effect and this translates in T (or X) being a function of everything: choice of boundary 
conditions, choice of action (both gauge and fermionic), volume, etc... The dependence of this 
parameter on these details is not small (see for example Fig. 3, where T changes by 100% when 
the renormalized coupling changes in the range ~ 2 — 4). In the end T (or X) cannot be computed 
analytically nor by studying other flow quantities, but has to be inferred from the same data that 
one is trying to improve. 

3.3 Symanzik improvement and the Zeuthen flow 

A very appealing approach to systematically reduce cutoff effects is the Symanzik improve¬ 
ment program. Its purpose is far more ambitious than looking for an accidental cancellation be¬ 
tween cutoff effects for some particular observable. When we choose a particular discretization (to 
define an acfion/observable) we are choosing what our theory looks like at very high energy scales 
(the cutoff I/a). The idea of the Symanzik improvement program is to choose these discretizations 
so that the effective theory at energy scales much smaller than the cutoff looks as close as possible 
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to the continuum theory. This program starts by describing a particular lattice action S'latt by an ef¬ 
fective continuum theory. For example for pure gauge theories any lattice action can be described 
by an effective action given by 

5eff = 5 ym + a^S2 + a^S4 + ..., (3.11) 


(i.e. any choice of actions will have the Yang-Mills action as first approximation S'ym. but different 
discretizations will produce different cutoff effects described by 52 , 54 ,...). Symanzik’s idea was 
to choose a discretization such that the effective continuum description has some of these terms 
absent (i.e. S 2 = 0 ). 

We can say that the basic idea behind Symanzik improvement program is to help us choosing a 
discretization that represents better continuum physics. Being a choice made at the cutoff scale I/a, 
it should be obvious that improvement parameters can not depend^ on the boundary conditions, the 
volume, or the choice of observable^. In other words, if one finds fhe Symanzik improved acfion, 
fhis is valid for all observables (fhis includes all values of fhe paramefer c fhaf defines fhe coupling), 
all volumes, all boundary condilions,elc... 

The applicafion of fhe Symanzik improvemenf program fo fhe gradienf flow has fo be done 
using a 5D local field fheory [ 6 ], in which fhe flow lime becomes fhe fiffh coordinafe wifh values 
in fhe posilive real axis. The acfion in fhe continuum Euclidean space can be written as 


^cont ^ _ 1 I -if dt f d\Tr{Lf,ix,t)[dtBf,{x,t)-DyGy^]} , (3.12) 

z^o J Jo J 

where Pgv and G^v are fhe field sfrengfh af f = 0 and t > 0 respecfively. L^(x,t) is a Lagrange 
multiplier in fhe sense fhaf fhe pafh infegral over fhis field imposes fhe flow equation ( 1 . 1 ) in fhe 
bulk of fhe 5D space (i.e. for t > 0). 

The main poinf is fhaf due fo fhe classical nafure of fhe fheory af t > 0, one can complefely 
eliminale all (i.e. fo all orders in fhe coupling) ff{a^) cutoff effecls in fhe bulk [27]. Making a long 
sfory shorf, fhis non-perfurbafive improvemenf is achieved by using fhe Zeuthen flow fo integrate 
fhe flow equations 


a^d,V^{xf) = -gl 


a 

1 + 


'^d:,SMV) 


Vg{x,t), 


(3.13) 


where and D* are fhe forward and backward lattice covaiianf derivatives, and 5 lw is the tree- 
level Luscher-Weisz improved action (i.e. Eq.(3.2) with cq = 5/3,ci = —1/12,C 2 = C 3 = 0). In the 
same way one can define clasically improved discretizations for observables (like ^(v,?)), fhaf will 
nol infroduce any ff{a^) discretization effecls. 

In summary, from fhe fhi'ee discrefizafion choices fhaf one has to make in order fo sfudy any 
flow quantify, fhe Symanzik improvemenf program gives a way fo fix fhe discrefizafion used fo 
infegrafe fhe flow equafion and fo evaluate flow observables, so fhaf no ff{a^) effecls are inlroduced. 
The remaining cutoff effecls are only fhose of fhe action used fo produce fhe configuralions and a 
couple of boundary counferferms. 

^Up to small logarithmic corrections. 

^One should nevertheless use an improved definition for the observable if one wants cutoff effects be reduced for a 
particular quantity. 
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Figure 4: Step scaling function in the minimal walking technicolor in two different schemes. Left: 
a2(gsp)/gsF scheme, showing an IR fixed point at g|p 2. Right: a 2 (g^p)/g^p in the 

gradient flow scheme, showing an IR fixed point at gQp ~ 5. (Source [34]) 


4. Some recent results 

Since the proposal of using the gradient flow to define a renormalized coupling there have been 
several works using this technique to investigate the dynamics of strongly coupled non-abelian 
gauge theories. Most of the applications are related to models beyond QCD and the search for an 
infrared fixed-point (see for a recent review [33]). On the other hand, the gradient flow has been 
used in QCD mainly for scale setting (see [8]). We have already mentioned some of the recent 
works in the previous sections. In this section we will briefly mention other works paying attention 
to the details discussed in earlier sections. The reader interested in a more exhaustive review of the 
recent works should also consult [33]. 

The work [34] studies the running of the coupling in the minimal technicolor model (SU (2) 
gauge group with two adjoint fermions) using both the Schrddinger functional coupling and the 
gradient flow coupling. Both schemes show an IR fixed point (see Fig 4) at quite different values of 
the coupling. Of course since the schemes are different there is no reason to expect the values of the 
coupling at the fixed point to be the same. This work uses a smeared gauge action, which modifies 
the lattice leading order value of t^{E{x,t)). The improved coupling definition Eq.(3.6) cannot be 
used without an analytic computation. The authors address the issue of large cutoff effects by using 
a large value of c 

Computing non-perturbative effects in SU{N) gauge theories in the large N ’t Hooft limit 
(keeping A = g^N constant) via lattice techniques is a difficult task because the number of degrees 
of freedom that have to be simulated on a computer increase with N. Nevertheless, using the ideas 
of volume independence (see [35] for a recent review), the authors of [29] have computed the 
running of the t’Hooft coupling Fig. 5. The authors use the tree-level corrected coupling definition 
Eq. (3.6) and the clover and plaquette definitions of the observable to estimate the systematic of the 
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Lattice Artefacts: Z(u}/u at u=0.98750 


Running coupling vs renormalization scale 
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Figure 5: Left: Continuum extrapolation of the step scaling function of the t’Hooft coupling in 
SU (oo). Right: Running of the t’Hooft coupling A = Ng^ in SU (oo) gauge theory (source [29]). 


continuum extrapolations. Their results also indicate that the tree-level corrected definition of the 
coupling Eq. (3.6) has a big impact in the cutoff effects. 

5. Fermion flow and renormalization of composite operators 

The flow can also be introduced for the fermion fields [6] 

, (5-1) 

where +B^, and the fermion flow field obeys the initial condition x{x,t)\t=o = 

this case the fermion flow field requires renormalization, but there is no operator mixing at positive 

flow time. i.e. composite fermion operators renormalize multiplicatively at positive flow time. 

In fact the study of Ward Identities at positive flow time has lead to a novel approach to the 
renormalization of some composite operators. We have recently seen the cases of the axial current, 
the chiral condensate [6, 36], or the energy momentum tensor [37]. 

On the other hand the so called small flow-time expansion offers the possibility of a general ap¬ 
proach to the renormalization of a generic composite operator 0(x,t). When the flow time is taken 
to zero 0(x,t) has to converge to the bare operator 0{x). Since the latter usually requires renor¬ 
malization, 0{x,t) becomes singular when taking the r — )• 0 limit in a way that can be expressed as 
an asymptotic expansion in renormalized operators with singular coefficients [5] 

0{x,t) = Y,Ca{t)0^{x) + ffit). (5.2) 

a 

The previous sum runs over all the operators of dimension smaller that the dimension of 0{x,t) and 
compatible with its continuum symmetries. We can use E{x,t) = —^Tr{G^vG^v), as an example, 
whose small flow time expansion reads 


E{x,t) = ci{t)l-\-C 2 {t){E^vF^v}R{x) + ff{t ). 


(5.3) 
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In this example, E{x,t) mixes with the vacuum and with the renormalized action density. In fact 
this relation can be used to determine the spin-0 component of the energy-momentum tensor [38] 

{T^g}R{x) = {F^^F^^}r{x) = limc 2 \t) {E{x,t) - {E{x,t))} . (5.4) 

This relation holds in the continuum, and its use in the lattice framework requires that the 
continuum limit is taken before the ? —)• 0 limit. This means that the scale given by the flow time has 
to be large compared with the lattice spacing, and small compared with any other relevant scale of 
the problem (i.e. a < 

Moreover one needs to compute the small flow time coefficients Ca{t). Ref. [38] proposes to 
use perturbation theory. This proposal has recently been used to extract thermodynamic quantities 
(see the contribution of this conference [7] for more details). Generally speaking there might be 
concerns about using leading order perturbation theory at energy scales below the lattice cutoff 
1/a, and this is the reason why we have seen different works focused on the non-perturbative 
determination of the small flow time coefficients. The authors of [39] propose to determine non- 
perturbatively the quantity Ya{t) = logCce(t) and then construct Ca{t) by integrating using 
perturbation theory as initial condition. The functions 7a (0 can be extracted from correlation 
functions. For example in the case mentioned before we have 

72 (?) = -2t^log(£'(x,t)£(v,5))connected + '^(0- (5-5) 

Although the idea is appealing, and has the advantage that it can be generally applicable to any 
composite operator with a complicated mixing pattern (including fermion composite operators), it 
remains to be seen how well it actually works in a concrete numerical example. 

We have also seen how the small flow-time expansion can be used to compute the strange 
content of the nucleon [40]. In this case we consider the connected part of the scalar strange 
density in a nucleon state 

0™^(x,t) = Nss{x,t)N^ — {NN^){ss{x,t)), (5.6) 

where N is a. nucleon interpolating operator. The small flow time expansion of this operator reads 

(9™'’(x,0 = C3(0O™‘’(x,0) + ff{t), (5.7) 

and c^{t) can be determined from the Ward identities (see [40] for more details). 

6. Conclusions 

We have seen in the recent years a growing interest in the gradient flow as a tool for scale 
setting and to define renormalized couplings. In this contribution we have reviewed the main ideas 
behind these applications putting a strong emphasis on the role of the continuum extrapolation 
and cutoff effects of flow quantities. We have also commented on the use of the small flow-time 
expansion as a general tool to renormalize composite operators, a very interesting and promising 
application where there is still a lot of work to do. 
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I would like to conclude with a brief comment on some recent works that I did not have 
time to cover or did not fit in the presentation, but nevertheless I found interesting. Perturbative 
computations have proven to be a useful tool (i.e. to understand the structure of cutoff effects), but 
also perturbative computations of flow quantities are usually tedious. The application of numerical 
stochastic perturbation theory for flow quantities [41] is an interesting idea. The use of the gradient 
flow to extract matrix elements non-perturbatively with a simplified mixing pattern has also been 
considered recently [42], and shows a promising future. 
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